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Elliptic flow of charm hadrons is investigated based on the quark coalescence model. Due to the 
large difference between the charm quark and light quark masses, hadrons containing both light 
and charm quarks show a qualitatively different V2{p±) from hadrons containing only light quarks. 
Simple relations are proposed to infer quark elliptic flow from those of hadrons. The effects of the 
finite momentum spread of hadron wavefunctions are also studied, and are found to be small for 
charm hadrons. 
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I. INTRODUCTION AND CONCLUSIONS 

During the early stage of ultra-relativistic heavy ion 
colhsions at the Relativistic Heavy Ion CoUider (RHIC) 
and the future Large Hadron Colhder (LHC), a decon- 
fined phase, often called the quark-gluon plasma, is ex- 
pected to be formed. As the dense matter expands, the 
parton degrees of freedom in the deconfined phase con- 
vert to hadrons. Quark coalescence is a simple model of 
hadronization, where it is assumed that gluons, which 
are abundant at high temperatures, either convert to 
quark-antiquark pairs or serve to "dress" quarks into con- 
stituent quarks near hadronization. The effective degrees 
of freedom are constituent quarks and anti-quarks, and 
mesons form from a quark-antiquark pair, while baryons 
from three quarks according to their valence quark com- 
positions. Quark coalescence has been applied to heavy 
ion collisions in the algebraic coalescence rehadroniza- 
tion (ALCOR) and microscopic coalescence rehadroniza- 
tion (MICOR) models Q to describe hadron abundances 
and in a multiphase transport model (AMPT) with string 
melting 2] to describe the elliptic flow at RHIC. Recently 
it has been used to address the elliptic flow of hadrons of 
different flavors and the large p/tt ratio [1110,13 

observed at RHIC. 

Elliptic flow, V2 = (cos(20)), the second Fourier mo- 
ment of the azimuthal momentum distribution, is one of 
the important experimental probes of collective dynam- 
ics in A + A reactions 0- It results from the spatial 
asymmetry in the transverse plane in non-central colli- 
sions, which is largest at early times. Therefore, elliptic 
flow is sensitive to the properties of dense matter, such 
as its equation of state .lO, 11, JJj or the effective 
scattering cross section of partons produced in the col- 
lisions Measurements of elliptic flow at 
high transverse momentum [l^ provide important con- 
straints about the density and effective energy loss of 
partons [11 d. 

Charm quark elliptic flow holds the key to resolve a 
recent controversy regarding charm dynamics. It was 
pointed out in Ref. jl9] that both i) initial perturbative 
QCD charm production without final state interactions 
(as modeled in PYTHIA 20]) and ii) complete thermal 
equilibrium for charm hadrons are consistent with the 



single electron spectra from PHENIX zf\. Information 
on the magnitude of charm quark elliptic flow is essential 
to distinguish between these two extreme scenarios. 

In this study we determine the elliptic flow of charm 
mesons and baryons based on quark coalescence and in- 
vestigate how to extract the charm quark elliptic flow 
from that of charm hadrons. We extend the simple quark 
coalescence formalism to the case of different constituent 
quark masses, which is relevant for heavy flavors, and 
also consider the effect of the finite momentum width of 
hadron wavefunctions. 

We find that relations between hadron and parton 
elliptic flow are modified for hadrons with constituent 
quarks of unequal masses. As a result, charm hadron 
elliptic flow is predicted to increase much slower with 
transverse momentum and saturate at a much higher mo- 
mentum scale than the flow of hadrons with only lighter 
flavors. A very slow increase of charm hadron elliptic 
flow with p± is a strong indication for a small charm 
quark elliptic flow. The effect of the finite momentum 
spread in hadron wavefunctions is relatively small for 
charm hadrons. In this case, simple relations exist to 
unfold both light and charm quark V2{pi_) from hadron 
elliptic flow. Comparing these predictions with upcoming 
experimental data will tell to what degree charm quarks 
rescatter in the medium and provide valuable insights 
into the dynamics and hadronization of the dense par- 
tonic matter in heavy ion collisions. 



II. MASS AND WAVEFUNCTION EFFECTS IN 
QUARK COALESCENCE 

A convenient starting point to describe meson produc- 
tion via coalescence af3 M is the relation [2^ 



E 
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dMlpp{q}\ fa{pa,x)fp{pi3,x) (1) 



between the phase space distributions of constituent 
quarks a and /3, the meson wavefunction ipj^ ^23*1 , and the 
invariant momentum distribution of produced mesons. 
Here p = jTa + p}), q = lhx— p h, pjvt} ^ the statistical fac- 
tor for the meson formation 7, 8, 22], and the integration 
runs over a 3D space-time hypersurface parameterized by 
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{x) . The expression is valid if coalescence is a rare pro- 
cess. When the coalescence probabilities are high for the 
constituents, instead of the quadratic (cubic) scaling of 
the meson (baryon) number with the quark number im- 
lied by the above equation, the scaling should be linear 
"H. Eq. CEl also neglects the hadronic binding energy. 
Because of the large charm quark mass, this assumption 
is better met for charm hadrons than for lighter ones. 
For quark coalescence into baryons a/37 ~^ Eq. 
can be generalized in a straightforward manner. 

Note, Eq. (Q) is valid for collectively expanding sources, 
provided the parton distribution functions, fi{p,x), 
change little over the spatial size of the hadron wave- 
functions (~ 1 fm). However, it has to be modified if flow 
velocity gradients are very large as shown by Eqs. (3.19- 
20) in Ref. . We leave the discussion of this and other 
space-momentum correlation effects to further studies. 

For a hadron at rest, has a small momentum space 
extension of Aqcd ~ 1/fm (based on the uncertainty 
principle) and therefore Eq. reduces to the simple for- 
mula considered in Ref. 5] . For a fast moving hadron, on 
the other hand, the wavefunction can change significantly 
due to Lorentz boost. We estimate this effect assuming 
that ipp is dominated by the contributions of dressed va- 
lence quarks. In this case, e.g., for mesons, |V'o('?'')P 
is the probability density for finding quark/ antiquark a 
and P with momenta (fc',?^') = (+9 V^: ^^V^) in the 
hadron rest frame (primed quantities refer to the hadron 
rest frame throughout this study). For a hadron at mid- 
rapidity with momentum p = {pj_nj_,0) in the LAB 
frame, the transverse momenta of valence quarks along 
the transverse boost direction n± are given by 



momentum fractions in the LAB frame are 



Szi = 



m 



H 



rriH 



(4) 



if the hadron is moving relativistically. For massive 
hadrons, such as charm hadrons, Szi ^ ql' ■ n±/mH ^ 1 
because ql ' is typically on the order of Aqcd • On 
the other hand, for hadrons with masses comparable to 
Aqcd, e.g., pions, Szi is always on the order of unity 
regardless of the transverse momentum. However, note 
that for such light mesons the binding energy cannot be 
neglected and hence the conventional coalescence formal- 
ism may be not be reliable. 

Since the constituent quark momentum components 
perpendicular to the hadron momentum (i.e., n± at mid- 
rapidity) are small (~ Aqqd), we further simplify Eq. 
by considering only the quark momentum component 
along n±. Then the integrals over the wavefunction can 
be recast, for mesons for example, as / d!^q\4'p{<l)\'^ = 
J dzal^piza)]"^ ■ As seen from Eq. |@}, Szi is independent 
of p_L for large p± and thus ^f{za) only depends on Zq. 
In this case Eq. Q becomes, with Zfj — 1 — z^, 



dN{p) fda^p^ 



d^p 



(2^) 



dZa\^M(Za)\ faiZaP,x)fi3(zpp,x). (5) 



III. ELLIPTIC FLOW FOR CHARM HADRONS 

A. Narrow wavefunction case (the limit of zero 
momentum spread) 



P±i 



A. 

uim' 



-P± + Qi ■n± 



ruM 



(2) 



Here E'^ = {mf + |gi'P)^^^, rm is the effective mass of 
valence quark i, and ttim is the meson mass. 

Suppose now that the rest-frame wavefunction was suf- 
ficiently narrow so that one can take q' 0. For the 
weakly bound system assumed, uim ~ rn-a + mp and 
Eq. ^ then gives p±i = p±mi/ (ma + mp). This re- 
lation also holds in general for the average constituent 
momenta, provided (g/)^ <C raf. Introducing the con- 
stituent momentum fractions Zi = p±i/p±, for mesons 
and baryons we then have, respectively, 
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(3) 



If the effective masses of constituent quarks are similar, 
P±a = P±/3 = P±/2 for a meson, which is the case, e.g., 
for pions or the J/V'- On the other hand, for D mesons 
TOq, <ti mfj and therefore most of the hadron momentum 
is carried by the heavy quark. The asymmetric momen- 
tum configuration arises because coalescence requires the 
constituents to have similar velocities, not momenta. 

The Lorentz boost also affects the width of the hadron 
wavefunction. From Eq. ^ the spreads of valence quark 



Because for charm hadrons the spread Szi is small, it 
is a good approximation to consider only the mean quark 
momentum fraction Za in Eq. |(SJ). For simplicity, we shall 
neglect the spatial variation of fi on the hypersurface 
and assume that in non-central heavy ion collisions the 
cos(20) component is the only non-trivial term in the 
quark azimuthal distribution, i.e., 

MpD = (27r)3^ = h,ip^) [1 + 2v2Ap±) cos(20)] . (6) 

Eq. © then relates V2 {p± ) of mesons or baryons to those 
of the constituent quarks as 



{PA 



V2,. 



- V2,f3 



1 + 2 V2,a V2,p 



' V2., 



B, ^ 1'2,Q + W2,/3 + l'2,7 + 3 U2,Q W2,/3 ^2.7 

V2{p±) = ■ — — 7- 

1 + 2 {V2,a V2,I3 + V2,a ^2,7 + f2,/3 ^2,7) 

— + Wj3 + 



(7) 



where V2J = W2,i(^iP_L) is the elliptic flow of valence quark 
i at its average momentum. As a result, for hadrons 
consisting of both light and charm quarks, elliptic flow 
at a given p±^ contains the light quark elliptic flow at a 
much smaller p_L . 
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If the light quark density is so high that the coalescence 
probability for a charm quark is unity, Eq. ((TJ breaks 
down, and charm hadrons just inherit the flow contribu- 
tion from charm quarks Preliminary experimental 
data indicate that p±_ scale for light quarks below which 
this occurs might be as low as 1 GeV Tr| . Also note that 
our study does not include the independent fragmenta- 
tion of partons, which eventually dominates over coales- 
cence at very large 0, 0. I3> ^'^^ power-law 
shape of the parton spectra at high p± . Since the frag- 
mentation function of heavy quarks in vacuum is much 
harder than that of light quarks, the p±^ scale above which 
independent fragmentation starts to dominate can be dif- 
ferent for charm quarks than those for light quarks. 
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FIG. 1: Comparison of v-2,(j>i) for quarks and charm hadrons. 

Fig. n shows the mass effects on elliptic flow of charm 
hadrons, considering cither t'2,c(p±) = f2,q(p±) (unla- 
beled) or V2,cij>\_) = 0. For illustration, dashed curves 
correspond to a light quark flow i'2,i}(p_l) which increases 
linearly then saturates when p\_ > po- Solid lines instead 
represent a more realistic scenario, where W2.g(p±) = 
v^""^ tanh(p_L/po), as found in the Molnar's parton trans- 
port model (MPC) For J/^ and protons, which 
consist of quarks with similar masses, vi{j)\_) saturates 
at ^ n^Po at a value nyV^"'^ 0, where n„ = 2 and 
3 for mesons and baryons, respectively. However, for D 
mesons and Ac baryons elliptic flow increases with pj_ 
very slowly due to the large difference in the constituent 
quark masses. The reason is that light quarks in these 
hadrons carry only a small fraction of the hadron momen- 
tum. The slow increase is most pronounced when charm 
quarks have no elliptic flow. Furthermore, when charm 
quarks have elliptic flow, V2{p±) of D mesons and A^, 
baryons both exhibit two different slopes because the flow 
of the heavier constituents enters the saturation regime 
much earlier. In both cases, V2{p±) of these hadrons 
saturates at a scale well above Uvpo- These qualitative 
features could be tested by experiments at RHIC and the 
LHC, and they will provide key information on the dy- 



namics of charm quarks /hadrons in dense matter created 
in heavy ion collisions ,2^ ,2^ . 



B. Unfolding quark elliptic flow 

When Eq. is valid and the effect of the momentum 
spread is small, it is possible to unfold the quark elliptic 
flow from hadron elliptic flow using Eq. Q . For example, 
if V2{p±) is known for D and A^, we have 

v'2iP±) = ((2 + rK) - ((1 + rK) , 



2vi, 



P± , (8) 



where r = mc/mq is the ratio of the charm quark and 
light {u and d) quark effective mass. 



C. Numerical estimates including the momentum 
spreads of wavefunctions 

For rough numerical estimates, we consider the condi- 
tions in semi-peripheral heavy ion collisions at the top 
RHIC energy, = 200^1 GeV. With the assumption of 
Eq. ijHJl, one obtains 

^~(Pp~ J dZa\^MiZa)f ha{ZaP±)hj3{zi3p±) 

X [1-H2w2,„(2„p_l)cos(20)] [H-2v2,/3(^/3PJ_)cos(2(/))] . (9) 

We recover Eq. when the momentum spread is zero, 
i.e., I'&A/l-Za)! = ^i^a — Za). Including the momentum 
spread but neglecting small corrections of higher orders 
in 112,4 7 the meson elliptic flow is 

M/ N SdZaWM{Za,Pl)[v2.a{ZaPL)+V2,D{zi3PAy\ 
J dZaWM(Za,P±) 

where the weight function for mesons is given by 

WM{Za,P±) = \^M{Za)\'^ ha{ZaP±)hp{zpp±_). (11) 

For baryons, the similar expression for w^(p_l) involves 
two integrals, over z„ and z^, where WBiza, Z0,p±) con- 
tains the product of three quark distributions. 

For the hadron wavefunctions, for simplicity we take 
the form used in the valon model 6]: 

\<^>MiZc.)f (X Z^z'^p, |$s(Za,2^)|'0CZ^z;^Z^, (12) 

where zp — 1 — Za for mesons and z-y — 1 — — 
for baryons. The normalization constants play no role in 
Eq. (|l()|l for the elliptic flow and have been omitted in 
the above. We take Aqcd = 0.2 GeV, and constituent 
quark masses m„ = = 0.3 GeV, nis — 0.5 GeV, 
rric = 1.5 GeV. Observing Eq. Q in the large p± limit, 
we take^Szi^ = n„ * (Aqcd/'71m)^/3. With given by 
Eq. Q, the exponents a and b are then determined. With 
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the convention that a is for the hghtest and b (or d for a 
baryon) is for the heaviest quark in a hadron, the values 
of (a, 6) are (0.25, 1.1) for kaons, (8.2,8.2) for 0, (4.9,28) 
for D, and (88, 88) for J/tp. In general, a and b increase 
for heavier mesons, reflecting their narrower spread in z. 
For pions the above method gives a singular solution, i.e., 
a = b < —1 leading to divergence in / dza\^Miza)\ . 
Therefore, for pions we take a = 6 = @, i.e., a flat 
distribution in z, for simplicity. For baryons, the (a, b, d) 
values are (3.6,3.6,3.6) for protons, (4.2,4.2,7.7) for A, 
(5.3, 9.5, 9.5) for S, (14, 14, 14) for ft, and (7.2, 7.2, 40) for 
A,. 
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FIG. 2: Weight function w{za,p±) of tt and at different p±. 

For heavy ion collisions at RHIC energies, the par- 
ton momentum distributions at mid-rapidity, hi{p_i), 
may be parameterized by the sum of an exponential 
soft part and a power-law hard component. Assuming 
these have identical shapes for light quarks (m, d, u, d, s 
and s), we take hi{p±) oc 6"^-^/^ + ch/{1 +p±/^Hy" 
with T = 0.18 GeV, ch = 0.36, A^ = 1.3 GeV, and 
w = 8.3. For charm (and anti-charm) quarks, we use 
hc{p±) oc (p_L-H0.5GeV)V(l+pj./6.8GeV)2i, which rea- 
sonably parameterizes the transverse momentum spec- 
trum of primary charm quarks at mid-rapidity in pp col- 
lisions at Vs = 200 GeV from PYTHIA 20]. We note 
that medium effects for charm quarks in A + A collisions 
such as energy loss |2^ |2^ are not considered in this 
study. 

Fig. 121 shows the weight functions l|ll(l for pions and 
(j) mesons at pj^ = 1, 4 and 8 GeV, with all maxima 
normalized to one. For the pion, they are not centered 
around the mean value Za = 0.5 at higher p±^ values. In- 
stead the dominant momentum configuration for the two 
constituent quarks is around Zq, ~ and 1. This asym- 
metric momentum configuration at high pi_ for wavefunc- 
tions that are flat in z is a result of the power-law en- 
hancement of the parton spectra at high p± , as shown in 
Ref. 0. For massive hadrons, such as the 0, the wave- 
function becomes narrower. Actually, for large a and 



b the meson wavefunction in Eq. 112|) becomes a Gaus- 
sian, exp[—{za — Za)'^ / {26z'^)], where Za — a/{a + b) and 
6z'^ ~ ab/{a + b)^. The narrowness of the wavefunc- 
tion then dominates over the momentum distributions in 
Eq. I|ll(l so that symmetric momentum configurations are 
usually favored. Indeed, as shown in Fig. |21 the weight 
functions for (f) center at Za = 0.5 for all three p± values. 
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FIG. 3: Elliptic flows of hadrons including momentum spread 
in wavefunctions: a) for mesons and b) for baryons. 

We now study the effect of momentum spread on 
hadron elliptic flow via Eq. (|10|l for mesons and the 
analogous equation for baryons. For light quarks in- 
cluding s and s, the elliptic flow is parameterized as 
W2,g(p_L) = tanh(p^/po) with = 0.05 and 

po = 0.75 GeV, based on parton transport simulations 
using MFC. We consider W2,c(pj_) = 'Vi,q{pi.) unless spec- 
ified otherwise. As shown in Fig. |3| all qualitative fea- 
tures for charm hadrons discussed earlier remain valid. 
However, in general hadron elliptic flow is reduced rela- 
tive to Eq. |7J) because the concave shape of the quark 

(j3±) ansatz tends to penalize any momentum spread 
(the average of vi (p_L ) at two different p j_ values is lower 
than the value at the average p_l). For example, the value 
of elliptic flow at = 6 GeV is lower than Eq. by 
36% for pions, 32% for kaons, 11%(24%) for D mesons 
with(without) charm flow, 9% for protons, 10% for A, 
5% for S, 10%(17%) for A^ with(without) charm flow. 
These corrections are especially large for pions and kaons 
at high p\_ and cause the decrease in their i'2(pj_) above 
~ 2 GeV seen in Fig. 13 The reason for this suppression 
is that the elliptic flow of the slower constituent quark of 
a pion or a kaon is still far below the saturation value be- 
cause of the very asymmetric momentum configuration 
(see Fig. However, note that the flows of pions 
and kaons may be significantly modified by resonance 
contributions and the binding energy, which have been 
neglected in this study. We also found that the momen- 
tum spread has a negligible effect on the elliptic flow of 
massive hadrons, especially those with quarks of simi- 
lar masses. The flow reductions are 0.8% for ^, 0.01% 
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